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Abstract: We consider antisymmetric perturbations of real symmetric matrices in the 
context of random matrix theory and two-color quantum chromodynamics. We investigate 
the level spacing distributions of eigenvalues that remain real or become complex conjugate 
pairs under the perturbation. We work out analytical surmises from small matrices and 
show that they describe the level spacings of large random matrices. As expected from 
symmetry arguments, these level spacings also apply to the overlap Dirac operator for 
two-color QCD with chemical potential. 
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1 Introduction 

Random matrix theory (RMT) has a remarkably wide range of applications in mathematics, 
physics and beyond. It describes universal statistical quantities that are not determined by 
the specific dynamics of a given system, but rather by its anti-unitary symmetries. One of 
the most prominent RMT quantities is the level spacing distribution P{s), the probability 
of two neighboring eigenvalues to appear at a distance s, which is universal after unfolding 
(to refer to an eigenvalue density normalized to one). For the classical Gaussian RMT 
ensembles, i.e. hermitian matrices with Gaussian probability distribution of the elements, 
the spacing distributions can be evaluated for any matrix size [1] . The results are very close 
to the level spacings of small matrices, a fact which is known as Wigner's surmise [2, p. 199]. 
The spacing distributions -P(s) of the smallest matrices are linear/quadratic/quartic in s 
times an exponential decay in (see e.g. Eq. (4.4)), where the power in the repulsion for 
small spacings - which equals the Dyson index (3 G {1,2,4} of these ensembles - counts 
the degrees of freedom per matrix entry in the ensemble. 
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In particular, RMT applies to quantum chromodynamics (QCD), where it describes 
the spectrum of the Dirac operator; for a review see [3]. The crucial feature is that QCD 
(with massless quarks) obeys a global chiral symmetry, which is spontaneously broken. In 
the e-regime, a finite volume regime dominated by the corresponding Goldstone bosons 
(pions), the equivalence with RMT has been proven in Refs. [4-6]. The appropriate RMT 
describing the microscopic spectral properties of QCD in this regime is the chiral ensemble. 
For bulk correlations, i.e. away from the origin, the chiral symmetry is no longer important 
[7-9] and universal properties of QCD spectra, such as level spacings, can be described by 
the standard Gaussian RMT ensembles [10, 11]. 

While the conventional RMT ensembles describe chaotic systems with definite sym- 
metries, physical systems are often 'not so ideal'. We briefly describe cases where the 
symmetry of the system deviates from this ideal situation. When no analytical expressions 
are available for the correlations in such transitions, one can have recourse to surmises. 

Firstly, physical systems may obey a certain anti-unitary symmetry only approximately 
or contain several parts with different anti- unitary symmetries. The latter leads to tran- 
sitions from one symmetry class to another. Including the Poissonian ensemble, which is 
the equivalent for generic integrable systems, the corresponding spacing distributions have 
been shown to fulfill a generalized Wigner surmise [12, 13]. Such transitions were observed 
at high temperatures in two- and three-color QCD [14-17] and are expected in the contin- 
uum limit of the staggered lattice Dirac operator in two-color and adjoint QCD [18, 19]. 
We will show a mixed-symmetry surmise at work for two-color QCD at imaginary chemical 
potential. 

Secondly, hermiticity can be broken in physical systems, as is the case for QCD with 
real chemical potential. In this case the determinant of the Dirac operator becomes com- 
plex, which hampers progress in numerical simulations of this system. Non-hermitian 
Gaussian random matrices were first discussed by Ginibre [20] . For the associated complex 
eigenvalues one can define a nearest neighbor individually and study the corresponding 
spacing distributions. For the complex Ginibre ensemble the level spacings do not fulfill a 
surmise from 2x2 matrices (even though they are always cubic for small s) [21, 22]. The 
spacings of this ensemble were compared to QCD spectra and showed agreement only for 
a small range of the chemical potential [23]. 

In QCD the chemical potential drives the breaking of the anti-hermiticity. An adequate 
random matrix model describing such a situation consists of two matrices where the break- 
ing of the hermiticity is governed by a real parameter [24, 25]. The eigenvalue densities and 
lowest eigenvalue distributions of such ensembles have been worked out in the microscopic, 
weakly non-hermitian limit (where the microscopic eigenvalues zN and perturbation iJ?N 
are kept fixed as the matrix size — t- do) [25-28] and successfully compared to quenched 
QCD data [28, 29]. 

The treatment of real random matrices is known to be an especially arduous task. The 
joint probability function of the eigenvalues for the real asymmetric ensemble (also called 
the real elliptic Ginibre ensemble) was calculated in [30] and generalized to the chiral case in 
[31]. Only recently, the correlation functions were worked out for the real Ginibre ensemble 
[32-34], the real elliptic Ginibre ensemble [35] and the real chiral ensemble [31, 36, 37]. 
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In this work we investigate the use of surmises to describe the level spacings of such real 
random matrix ensembles. More specifically, we examine the case where the unperturbed 
matrices are real and symmetric, i.e., they belong to the Gaussian orthogonal ensemble 
(GOE, /3 = 1), and the perturbations are real and antisymmetric. The perturbations 
are assumed to be small, such that the imaginary part of typical eigenvalues is at most 
of the order of the mean separation between neighboring eigenvalues along the real axis 
(whereas in the Ginibre ensemble hermiticity is maximally broken). For real matrices 
the flow of eigenvalues is somewhat special, because two real eigenvalues need to coalesce 
to 'give birth' to a complex conjugate pair, see below (this is not so for the complex 
(GUE, /3 = 2) or symplectic (GSE, /? = 4) ensembles: in the former the eigenvalues are 
not restricted to form complex conjugate pairs, but can be distributed arbitrarily in the 
complex plane; in the latter the unperturbed matrices have two-fold degenerate eigenvalues 
due to Kramers' degeneracy and these split up and become complex for arbitrary small 
anti-hermitian perturbation). Actually, the original GOE eigenvalues are attracted by an 
antisymmetric perturbation. For eigenvalues that meet in such a process and become a 
complex conjugate pair, a spacing can be uniquely defined as the distance between them. 
As we will show, the distribution of these spacings for large random matrices can be 
approximated very well by that of small matrices. The complex eigenvalues leave a trace 
in the statistics of the real eigenvalues, too. We will work out surmises for the spacings 
among the latter and show that they also apply to large matrices. We interpret these 
findings in the same way as for other Wigner surmises, namely that the dynamics of a 
neighboring pair is not much influenced by the effect of all other eigenvalues. 

As an application we analyze spectra of two-color QCD, i.e., with gauge group SU(2). 
As the latter is pseudo-real it ensues that the massless continuum Dirac operator for 
fermions in the fundamental representation belongs to the GOE universality class (see 
Sec. 6.1). This property holds for vanishing chemical potential. Although the chemical 
potential perturbs the anti-hermiticity of the Dirac operator, it keeps its anti-unitary sym- 
metry. Therefore, the random matrix results discussed above apply to two-color QCD 
at finite chemical potential, both in the continuum and on the lattice with the overlap 
operator. 

This paper is organized as follows. After describing the symmetries and corresponding 
spectra of matrices in the next section, Sec. 3 is devoted to the eigenvalue dynamics in 
perturbation theory. In Sec. 4 we work out the surmises for the different types of spacings. 
These are then matched to data from large random matrices with mixed hermiticity in Sec. 5 
and to lattice two-color QCD data with chemical potential in Sec. 6. We also measured 
spacings at imaginary chemical potentials and discuss the comparison to a mixed-symmetry 
surmise in the last subsection of Sec. 6. We summarize our findings and conclude in Sec. 7. 

2 The setting: symmetries and spectra of RMT ensembles 

Our unperturbed random matrices Hq shall be real and symmetric. We therefore choose 
them from the GOE, so that the entries are Gaussian random numbers with mean zero 
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Figure 1. Evolution of eigenvalues of a typical random matrix with increasing non-hermitian 
part (schematically). The top panel depicts a part of the spectrum for the unperturbed case, i.e., 
for a real symmetric (GOE) matrix. Towards the bottom panel the antisymmetric part increases, 
which results in some eigenvalues getting very close and finally becoming complex (the rightmost 
eigenvalues are spectators). Different types of level spacings are sketched: (1) spacings between 
original eigenvalues and those that remain on-axis and have no interspaced complex conjugate pair, 
named P^^^ and described by the pure GOE surmise, (2) spacings between on-axis eigenvalues with 
interspaced complex conjugate pairs, named Pq„, (off) spacings between complex conjugate pairs 
called off-axis, named Poff- Note that for QCD spectra the real and imaginary axes are swapped. 



and variance 

{{Ho)l) = l and {{Ho)l) = l ^ ^'^ ' ^^"^^ 

Gaussian random matrix ensembles can be characterized by whether they commute 
with an anti-unitary V = UK, where U is unitary and K is the complex conjugation 
operator 

[Ho,V] = equiv. to = U^HqU, (2.2) 

and whether V squares to 1 or —1. For the GOE one has 

= 1 equiv. to UU* = 1 , (2.3) 

from which it follows that the matrices Hq can be chosen real. 

In this paper we will study random matrices from the GOE undergoing a real anti- 
symmetric perturbation with a tunable coupling parameter. The characteristic equation of 
such a mixed matrix has real coefficients and consequently the eigenvalues are either real or 
come in complex conjugate pairs. Typical spectra of such matrices are depicted schemati- 
cally in Fig. 1 together with different types of spacings between them. In the following, the 
real eigenvalues will be called 'on-axis' and those with a non-zero imaginary part 'off-axis'. 
This notation also applies to the case of QCD, where due to the anti-hermiticity of the 
Dirac operator at vanishing chemical potential the axes are swapped (then on-axis refers 
to purely imaginary eigenvalues and off-axis to those evolving from there into the complex 
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plane) . Obviously the spacings between on-axis eigenvalues and off-axis eigenvalues will be 
treated separately. Among the former we will further distinguish between those with and 
without an interspaced complex conjugate eigenvalue pair. This distinction only makes 
sense if the hermiticity of the initial GOE matrix is only mildly broken (since otherwise 
the off-axis eigenvalues move arbitrarily in the complex plane and cannot be assigned to 
a single on-axis pair anymore), and we will therefore only consider this case here. This is 
also why we do not consider the rare case of two or more interspaced eigenvalue pairs. 

All our level spacing distributions are subject to the normalizations ds P{s) = 1 
and (s) = ds s P{s) = 1. 

3 Repulsion and attraction of eigenvalues in perturbation theory 

To gain intuition about the dynamics of the eigenvalues, we consider a real symmetric 
matrix Hq, with real eigenvalues 9i and corresponding real eigenvectors perturbed by 
a real matrix AM with A G M.q . We compare the effect of symmetric and antisymmetric 
perturbations M, whose matrix elements in the unperturbed basis, 

Mij = {^i\M\^,) , (3.1) 

are real and symmetric or antisymmetric in too, 

Mji = ±Mij , for = ±M . (3.2) 

As is well known, ordinary perturbation theory up to second order yields for the eigenvalues, 

9,^9, + XMu + A^ V . (3.3) 

The first order term induces a random walk of the eigenvalues without correlating them. 
Moreover, in the antisymmetric case such expectation values are zero, so there is no first 
order contribution. 

When omitting the random walk term, the eigenvalue differences are 



(9,-9j)^i9i-9,) 



M2. 



)2 yi^Oi - 9,){9i,, - 9k^,j) 



(3.4) 



If 9i and 9j are much closer to each other than to any other eigenvalue, the third term in 
square brackets can be neglected in relation to the second one. Then, for real symmetric M 
(upper sign) the difference grows, which is the common knowledge that eigenvalues repel. 

In contrast, for real antisymmetric M (lower sign), which we will consider from now on, 
the eigenvalues attract. Thus in perturbation theory it is favored that eigenvalues coalesce. 

We recall that our investigations so far are based on the eigenvalues 9i of the unper- 
turbed matrix, which are real and remain so as long as ordinary perturbation theory is 
valid. ^ 



^Ordinary perturbation theory only yields real shifts of the eigenvalues to all orders. 
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When two eigenvalues Oi and 6j are very close, one has to apply almost degenerate 
perturbation theory [38] . Neglecting the second order contribution of all other eigenvalues, 
it gives for the perturbed eigenvalues 



evs 



Oi \ I M- 



+ A 



Q 9j \ -Mij 



(3.5) 



where we have made use of Eq. (3.2), with Mij still given by Eq. (3.1), and evs denotes 
the eigenvalues of the matrix. For the eigenvalue difference one obtains 

s = ^{e, - - Axmf. . (3.6) 

To second order in A the expansion of S reproduces the result from second-order ordinary 
perturbation theory discussed above. 

More interestingly, the argument of the square root can turn negative for some critical 
A, above which the eigenvalue difference becomes purely imaginary. The latter is consistent 
with the fact that the complex eigenvalues of a real matrix come in complex conjugate pairs. 
Imaginary spacings occur for A bigger than 

(specific to each individual eigenvalue pair). In the exceptional case that Mij exactly van- 
ishes, this particular eigenvalue difference remains unperturbed for any A (in perturbation 
theory). For perturbations beyond Acrit the eigenvalue difference can be written as 



S = 2i\Mij\^\^ - Xl,, = 2i\M,jy5\{2\,,,t + 5\) , ^ = A - Acnt , (3.8) 

and increases with the excess 6\ of the coupling over the critical coupling. Therefore, the 
perturbation causes the eigenvalues to be repelled in the complex plane. 

Moreover, the eigenvalues leave the real axis perpendicularly to it. This can be seen 
by an explicit calculation of Oij in (3.5) for A > Acrit, but it also follows directly from 
the fact that the antisymmetric perturbation is traceless and thus does not change the 
center of mass of the eigenvalues. When considering antisymmetric perturbations of large 
matrices beyond perturbation theory, only the center of mass of all eigenvalues is fixed 
and a complex conjugate pair could in principle leave the axis at any angle, as its center 
of mass contribution could be compensated by the motion of the remaining eigenvalues. 
Our data confirmed that, in accordance with almost degenerate perturbation theory, the 
complex eigenvalue pairs move perpendicularly to the axis just after their creation. When 
the perturbation grows further, higher order effects set in and the complex conjugate 
eigenvalues can leave the trajectory perpendicular to the axis. 

4 Surmises 

The idea of Wigner surmises is to approximate level spacing distributions of large random 
matrices by those of random matrices with smallest possible size, typically 2x2 matrices. 
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Note that in the perturbative treatment in the previous section we have focussed on a 
pair of nearest neighbor eigenvalues and thus ended up discussing 2x2 matrices and the 
properties of their eigenvalues. This may be seen as an indication that surmises for mixed 
systems may apply to large random matrices. 

We recall that we consider three different types of spacings in our (large random matrix 
and QCD) systems: on-axis spacings with one and without interspaced complex conjugate 
eigenvalue pair, and off-axis spacings (see Fig. 1). For the first and third type we will 
deduce surmises from 2x2 matrices, computing the corresponding spacings. The spacings 
of the second type will be computed from the joined probability density for 4x4 matrices. 

4.1 Prelude: 2x2 matrices and their spacings 

For size 2 x 2 let us consider a real symmetric and traceless matrix perturbed by an 
antisymmetric one, 



^= 7 +A ^ : , (4.1) 




with random numbers a, b and c that are Gaussian distributed with zero mean and unit 
variance. Note that the matrix H with A = is equivalent to a 2 x 2 GOE matrix, up 
to a common shift of the eigenvalues which does not alter the spacings. In the GOE the 
diagonal entries have twice the variance of the off-diagonal ones, but it can easily be shown 
that neglecting the part proportional to the identity, as we have done in H to simplify the 
computations, renders the variances equal. The eigenvalues of H are ib\/a^ + b'^ — c^X^, 
which, for a particular draw of random numbers, turn imaginary beyond a critical coupling 

Aen, = . (4.2) 

\c\ 

For A < Acrit the spacing is 2V a'^ + b'^ — c^A^ and hence real, while for A > Acrit the spacing 
is 2i\/ c^A^ — a? — h"^ and hence purely imaginary. 

The probability that the matrix H has an imaginary spacing can be easily computed 
by integrating over all random matrix entries with the constraint that A is above its critical 
value, 

/oo -| 
dadhdce-^^'+^^+''''^/^e{c^\^ -a'-b'') = l- ^ . (4.3) 
-oo Vl + A^ 

As expected, the limiting cases yield p{X — )• 0) = and p{\ — )• oo) = 1. 

4.2 On-axis without interspaced complex eigenvalues 

For on-axis spacings without interspaced complex eigenvalues our surmise is 

PL{s) = lse-^^^\ (4.4) 

which means that even in the presence of a (small) antisymmetric perturbation these 
spacings follow the GOE surmise governing the unperturbed spacings. This is so because 
Eq. (4.4) is the distribution of real spacings^ of H independently of the coupling parameter 



^The plural refers to the ensemble, as one representative of this random matrix has only one spacing, 
which will only be included here if real. 
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A. Let us demonstrate the derivation of this distribution exphcitly, as a typical example. 
We start with the distribution of the non-normalized spacing S 



-oo 



2tt 



dc 

oo JQ 



dr r e-('''+^')/2 5 (s - 2\/r'^-c^\A 9 (r^ - c^A^) 

In the integration over r, the J-function replaces 

r ^ ro = V52/4 + c2A2 

with Jacobian 



d 
dr 



5-2x/r2-c2A2 



r=ro 



5 



5 



4v^52/4 + c2A2 4ro 



Therefore, we obtain 



Q{S) 



IT 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



The integration over c just gives a constant factor, and the ^-function forces the spacing S 
to be real as assumed. After normalization, we obtain P^^ from (4.4). 

4.3 On-axis with interspaced complex eigenvalues 

When two on-axis eigenvalues have an interspaced complex eigenvalue pair, we expect the 
latter to have a noticeable influence on the spacing of the on-axis eigenvalues. To obtain a 
surmise for the distribution of these spacings, we consider a 4 x 4 matrix 

H = Ho + XA, (4.9) 

with Hq taken from the GOE and an antisymmetric matrix A with probability density 

u;(^)~e-M^^^), (4.10) 

i.e. Gaussian distributed entries having the same variance as the off-diagonal entries of 
Hq. This is the smallest possible matrix that can yield a surmise for this case as there are 
four eigenvalues which are relevant for the dynamics of the spacing. The joint probability 
distribution of the eigenvalues of H depends on the number of real and complex eigenvalues. 
For two real eigenvalues 6i and 02 and one complex conjugate pair z and z* , it reads up 
to a normalization [36] 

P{9i,92, z, z*) oc e-^?-^i-^'-(^*)' erfc ^ ^^^"^^ - z*\^ 2iA{9i,92, z, z*) , (4.11) 

with the Vandermonde determinant 

A{9i,92, z, z*) = {z* - z){z* - 92){z* - 9i){z - 92){z - 9i){92 - 9i) , (4.12) 
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Figure 2. Spectrum of the 4x4 matrix used for the surmise for 



where it is assumed that 02 > Oi and Dim z > 0. To obtain the distribution of the spacing 
5 between 6i and 02, we set 62 = 0i + S and introduce new variables a = Die z — 9i £ [0, S] 
and b = 3m z £ [0, 00), cf. Fig. 2. This results in 



P{ei,9i + S,ei + a + ib,ei+a- ih) 



(4.13) 



(X S b 



\a - Sf + b'] {a' + b') e~'l~io,+s?-2ia+e.f+2i.- ^^^^ ( ^V^^^ 



We integrate out ^i, a, b and perform an irrelevant rescaling S — )• 25 to obtain the 
normalized spacing distribution 



non- 



/■oo f'lS f'co 

/ d5 / da/ deiP{0i,ei + 2S,ei + a + ib,9i+a-ib) 

Jo Jo J-00 

{V^e^%rf(5) [aVi + A2 {3X^ + 85^) + 4 (^1 + A2 - l) (4^4 _ g^2 ^ 3^ 
-85 + A2 - 1) (252 - 1) - aVi + A2] } = Q(5; A) . 



DC Se~ 



Properly normalized, the spacing distribution reads 

P^,{s;X) = CDQ{Ds;X), 
with normalization and scaling factors 

12^/3 



C 



and 



,00 \ -1 

J dSQ{S;X) 



XWl + X^ (9A2 + 8) + 8 (Vl + A2 - 1 



(4.14) 



(4.15) 



(4.16) 



00 

D = C J dS SQ{S;X) = ^ 26 (\/l + A2 - l) - 27\/2 arccot(V2) 


+Vl + A2 (^18A^ + 20A2 + 27\/2 (l + A2)^ arccot(\/2) 
In the limit A — t- 0, the spacing distribution reads 



(4.17) 



limp2^(s;A) = ^se 
A^-O 87r 



Ak^s^ - 6ks + V7i^e'''^'(4KS^ - 4^252 + 3) erf(Ks) , (4.18) 
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where 



The opposite hmit A — )• oo yields 



lim Pol(s; A) = 4\/3ry2se-2(''^) erf(77s) , (4.20) 

A— >oo 



with 



2 + 3\/2arccotV2 

„ = hm L> = ^ — w 0.7510 . (4.21 

A^oo 2^/3^ 

For small spacings s and non-zero A, the spacing distribution P^^ is proportional to 
s^, whereas in the limit A — >• the first term in the Taylor expansion is proportional to s^. 
One can understand these repulsion strengths from the Vandermonde determinant in the 
joint probability density alone, i.e. by focussing on the linear repulsion of every eigenvalue 
pair and neglecting the exponential and erfc factors (which only contribute significantly at 
larger spacings). The relevant integral is 



s 

da{z* - e2){z* - 9i){z - e2){z - ei){92 - Oi) 
fS 



z=0i+a+ib,02-ei=S 



= S [ da [{a- Sf + b^] {a" + h") = h^S^ + I b^S^ + ^ ■ (4-22) 
JO 3 30 

This is indeed proportional to for small S, unless the distribution of b becomes a 6- 
function around zero, which happens for A — )■ 0. Then, the distribution is proportional to 

Spacings for various parameters A are plotted in the left panel of Fig. 3. 
4.4 Off-axis 

To get a surmise for the spacings between complex conjugate eigenvalue pairs, we consider 
again the 2x2 matrix H given in Eq. (4.1). This time, we are interested in the distribution 
of the imaginary spacings of this matrix. In analogy to Sec. 4.2 we define Q{S) as in 
Eq. (4.5) with appropriate arguments of the delta and step function, 

/oo 
da db dc e-{"'+'''+^')/2 5(S- l^c^X^ - - bA 6 (c^A^ - a^ - b^) . (4.23) 
-co ^ ' 

For the normalized spacing distribution we obtain with calculations similar to those of 
Sec. 4.2, 

Poff(s; A) = CD^se-^'^'erfc {Osy/l + X^/)^ , (4.24) 

with constants 

C= ^ , = ^. (4.25) 
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Figure 3. Left: Wigner surmise P^^(s;X), Eq. (4.15), for on-axis eigenvalues with an interspaced 
complex conjugate eigenvalue pair for A = 0, 0.3, 0.5, 1, oo (decreasing maxima). Right: Wigner sur- 
mise Poff(s; A), Eq. (4.24), for off-axis spacings and couplings A = 0, 2, 10, oo (decreasing maxima). 
Plotted for comparison is the Wigner surmise for the GOE (dashed), which is the surmise for the 
on-axis spacings with no interspaced complex pair, Eq. (4.4). 

For small s the level spacing is linear, just as in the GOE. This is understandable, as 
the derivation of Poff involves an integral over one off-diagonal random number which is 
similar to that encountered in the derivation of the GOE spacing. The number of degrees 
of freedom, which governs the level spacing at small s, is thus the same for both cases. 
The limit for A — t- is 

64 / 4s \ 

hm Poff(s; A) = — s erfc . (4.26) 

Note that in this limit the perturbation is switched off. For statistical reasons, off-axis 
eigenvalues will exist for arbitrary small perturbations (see Eq. (4.2)). Equation (4.26) 
describes the (normalized) level spacings of these eigenvalues in the A — >• limit (whereas 
Eq. (4.18) reflects the influence of these eigenvalues on the neighboring on-axis spacings in 
this limit). Even though the difference between Eq. (4.26) and the GOE spacing distribu- 
tion is rather small, the two are clearly distinguishable in Fig. 3 (right). 
For A — 7- oo the distribution Poff is simply half a Gaussian, 

lim Pofr(s; A) = - exp ( --) , (4.27) 

A->oo vr \ / 

the spacing distribution of the perturbation alone. This limit, however, is not uniform 
at s = 0, since Pofr(0; A) = for every finite A, whereas lims_j.o limA-s.00 Poff(s; A) = 2/-K. 
Spacings for various couplings A are plotted in Fig. 3 (right) and this discontinuity is clearly 
visible. A similar effect has been observed for mixed symmetry classes of (small and large) 
random matrices in Ref. [13]. There, even a Gibbs-like overshoot of the curves near s = 
was observed, which is absent here. 
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We have checked that the surmises and Poff could equally well be obtained from 
the joint eigenvalue probability distribution of real 2x2 matrices [36], i.e. analogues of 
Eq. (4.11), with either two real or a pair of complex conjugate eigenvalues. 

With these surmises at hand, our aim will now be to describe the level spacing distri- 
butions of large random matrices, matching the parameter A to the corresponding coupling 
parameter and to QCD spectra, where the matching will be to the chemical potential. 



5 Comparison of the surmises to large RMT spectra 

To check the validity of the surmises calculated in the previous section, we applied them 
to the spectra of large dimensional random matrices of the form 

where Ho is real symmetric and taken from the GOE with matrix size A^, whereas A is real 
antisymmetric with probability density 

w{A) ~ e-'M^^n , (5.2) 

i.e. all elements of A are independently Gaussian distributed with the same variance as the 
off-diagonal entries of Hq (as was also used for the small matrices in the surmises, cf. e.g. 
Eq. (4.10)). The coupling parameter A is divided by y^lN/n in order to make it compara- 
ble to the one used in the surmises, A. As far as the spacing distribution is concerned, this 
normalization makes A a universal, A^-independent coupling parameter, see [13] for a de- 
tailed discussion. As also argued there, a constant spectral density is necessary to compare 
spacing distributions of large matrices from mixed universality classes to the correspond- 
ing distributions of small matrices. Otherwise, different coupling strengths are mixed. 
Therefore, we only evaluated eigenvalues with real part in the interval {—^/N /4), 
i.e. around the center of the real spectrum of H. As the spectral density is not exactly 
constant, we measured the on-axis spacings in units of the local mean spacing (obtained 
by an ensemble average), which is equivalent to unfolding the spectrum. For the off-axis 
spacings, no unfolding was done. 

The numerically obtained spacing distributions of H are shown in Fig. 4 for various 
values of the coupling parameter A. As can be seen, the surmises for the on-axis spacings of 
both types describe the data very well for coupling parameters up to A = 2. The coupling 
parameters A for the on-axis spacings of type 2 were obtained by a fit with least square 
minimization. 

For the off-axis spacings we are able to predict a 2 x 2 coupling parameter A through 
the frequency p of imaginary spacings: we measure the latter for the large matrices and 
choose A such that the same p is realized for 2 x 2-matrices, i.e. we invert Eq. (4.3) to obtain 
A = ^1/(1 -pf - 1. The surmises for the off-axis spacings obtained this way match the 
numerical data very well for A = 0.2 and 0.5, differ slightly for A = 1 and are far off for 
A = 2. Although the coupling parameter could also be determined by a fit to the level 
spacing data, we observed that this did not yield an improved estimate. Note that the 
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Figure 4. Spacing distributions of 400 x 400 random matrices of the form of Eq. (5.1), with 
various values of the couphng parameter A. Top: on-axis spacings of type 1, surmise given by 
GOE, Eq. (4.4). Middle: on-axis spacings of type 2, surmise given by Eq. (4.15). Bottom: off-axis 
spacings, surmise given by Eq. (4.24); the dashed curve is the GOE spacing for comparison. For 
each value of A, 2 • 10^ random matrices were diagonalized. 



surmises' maxima are always left of the maximum of the GOE, cf. Fig. 3 right, in contrast 
to the A = 2 maximum. Not surprisingly, the surmise for the off-axis spacing does not 
work for large coupling parameters, i.e. for large perturbations. The reason for the earlier 
break down of the off-axis surmise could be related to the additional degree of freedom of 
the eigenvalues when they move in the complex plane. 

6 Application to two color QCD with chemical potential 
6.1 Symmetry and hermiticity of the Dirac operator 

As far as anti-unitary symmetries are concerned, the surmises discussed in the previous 
sections should approximate the spacing distributions for two-color QCD with a real quark 
chemical potential, as will be analyzed now. 

6.1.1 Continuum 

In the presence of a chemical potential, the Euclidean-space Dirac operator in continuum 
two-color QCD is given by 

D = -fyDy + m + n'yi, (6.1) 
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with covariant derivative 



= du + iAl{x)Ta , (6.2) 

where are the Pauh matrices, the generators of SU(2). The massless Dirac operator at 
zero chemical potential is an anti-hermitian operator, which has purely imaginary eigen- 
values. The mass merely shifts the spectrum, while the chemical potential distorts it 
nontrivially. Since the former is irrelevant for level spacings, we will set m = for the 
following discussion. 

The anti-unitary symmetry of the SU(2) Dirac operator is based on the pseudo-real 
nature of SU(2) implying 

(iTaY = 'T2{iTa)T2 , (6.3) 

as well as on the charge conjugation properties of gamma matrices 

{il.T = C\i^,)C with C = 7274, (6.4) 
where we use the Weyl representation 

y{-iTi, I2) J 

Making contact to the real random matrices considered so far is easiest after multipli- 
cation with the imaginary unit, 

iD = i-fu (S) {dul2 + Al{x)iTa) + i/i74 (S) I2 , (6.6) 

where we have separated the spin and gauge parts explicitly. From Eq. (6.3) and Eq. (6.4) 
it follows that iD obeys the GOE symmetry 

{iD)* = U\iD)U , U = C0T2 with UU* = 1 . (6.7) 

As a consequence, iD can be made real by a unitary transformation. We denote such an 
equivalence by ~ and write 

iD^Ds + Da, (6.8) 

where we have split the operator into its real symmetric part Dg and real antisymmetric 
part Da- From the hermiticity at = one identifies iD{^. = 0) ~ Ds and further 
i/i74 ~ Da- Thus up to the factor i - which only rotates eigenvalues, but does not influence 
spacings - the situation for the continuum Dirac operator is the following: for vanishing 
chemical potential it is real symmetric with on-axis eigenvalues obeying the GOE symmetry, 
whereas the chemical potential introduces a real antisymmetric perturbation with possible 
creation of complex conjugate eigenvalues. This agrees with the setting for the random 
matrices with /i playing the role of the matrix coupling parameters A and A. 

The massless Dirac operator also satisfies chiral symmetry, {L',75} = 0, resulting in 
all nonzero eigenvalues coming in pairs ^iO. In the numerical results below we will restrict 
ourselves to the half of eigenvalues with positive imaginary part and consider spacings in 
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the bulk, away from 6 = 0. For the spacing distribution in the bulk the chiral symmetry is 
known to be irrelevant [7-11].'^ 

6.1.2 Lattice 

When simulating the theory on a space-time lattice, the Dirac operator has to be dis- 
cretized accordingly. Lattice discretizations of the Dirac operator typically alter some of 
its symmetry properties, such that its anti-unitarity properties can differ from those in the 
continuum^ and have to be analyzed carefully. 

To preserve the chiral symmetry of the discretized Dirac operator one uses the overlap 
fermion formulation [39-41] , whose definition, in the presence of a quark chemical potential 
H, is given by [29] 

Dov(^) = 1 + 75 sign(75Dw(/u)) , (6.9) 

where we have put the lattice spacing to unity, a = 1, sign is the matrix sign function 
satisfying (sign^)^ = 1, and 

3 

Z)w(/i) = 1 - K J] (7;+ + T;-) - K (e^T+ + e-^T^) (6.10) 

i=l 

with 

{T^)yx = (1 ± lu) Ux,±u Sy,x±u (6.11) 

is the Wilson Dirac operator at non-zero chemical potential [42] with hopping parameter 
K = l/(8+2mw), Wilson mass m„ S (—2, 0) and gauge configurations Ux,±u £ SU(2), where 
Ux-u = The exponential factors e^^ implement the quark chemical potential on 

the lattice. For /i = the argument of the sign function in Eq. (6.9) is hermitian, while for 
/i 7^ it is non-hermitian. 

To actually compute the overlap operator we need to define the matrix sign function 
for a general complex matrix A of dimension n. A generic matrix function f{A) can be 
defined by 

f{A) = ^.i f{z){zl - A)-Uz , (6.12) 
2711 Jy 

where F is a collection of contours in C such that / is analytic inside and on F and such 
that F encloses the spectrum of A. If A is diagonalizable, i.e., A = UAU~^, with diagonal 

^Note that chiral symmetry imphes another anti-unitarity relation, D* = {-ysUy D(j5U), which at first 
sight would mean that D belongs to the class of real antisymmetric operators. This is, however, not the case, 
since for the classification one first has to reduce the operator to irreducible blocks. This is best done on , 
in which chiral eigenvalue pairs diid become twice degenerate eigenvalues —8^. As commutes with 75, its 
eigenmodes can be chosen to be of definite chirality. The two relations (D^)* = D^U = (75 (7)^-0^ (75?/) 
coincide on the subspaces of fixed chirality, where both are proportional to w = r2 (8> r2 with mm* = 1. This 
qualifies and thus also D for GOE. 

''as is the case for the staggered Dirac operator [18, 19] 
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eigenvalue matrix A = diag(Ai, . . . , A^) and U G G/(n, C), then this general definition can 
be simplified to the well-known spectral form 

f{A) = C//(A)C/-i , (6.13) 

with 

/(A) = diag(/(Ai),...,/(A„)) . (6.14) 

If A cannot be diagonalized, a spectral definition of f{A) can still be derived using the 
Jordan decomposition [43]. For hermitian A the eigenvalues are real and their sign is 
defined by sign(x) = ±1 for x ^ with x E M, such that Eq. (6.13) readily defines 
the matrix sign function. For non-hermitian A the eigenvalues are complex and require 
a definition of sign(2:) for z G C The sign function needs to satisfy (signz)^ = 1 and 
reproduce the usual sign(a;) for real x. We define [44] 

sign(2;) = = sign(J)\ez) , (6.15) 

where the cut of the square root is chosen along the negative real axis. This choice, although 
not unique, gives the correct physical result for the overlap Dirac operator in Eq. (6.9) (see 
Ref. [44]). 

To investigate the anti-unitarity properties of the overlap operator we first observe 
that links Ux,±u as elements of the SU(2) group obey 

Ul^, = 4Ux,±,T2 (6.16) 

(as they can be written in terms of I2 and ir^). For the Wilson Dirac operator multiplied 
by 75 one can show that 

(75-Dvv)* = U^Dyj'j^)!/ (note the inverted ordering) (6.17) 

with U as in the continuum, see Eq. (6.7). 

For the overlap operator we have to defer the discussion of its anti-unitary symmetries, 
because it has no definite hermiticity, even at zero chemical potential. Instead it fulfills 
the Ginsparg- Wilson relation [45], 

{^ov,75} =^ov75^ov, (6.18) 

which, together with the 75-hermiticity, 75D0V = -Dov75, valid at vanishing /i, forces the 
eigenvalues on a circle 9 = 1 + exp(icp) in the complex plane. When investigating the 
spectrum of the overlap operator, we will rather consider a related operator [46] 



2I)ov _ 2 1 + 75 sign(75Dw) 
2 - Dov 1 - 75 sign(75Dw) ' 



= -^iig^ = 2 : ^ :;!:;.^. ; , (6.19) 



which satisfies exact chiral symmetry, {-Dp, 75} = 0, at the price of being non-local. For 
/X = this operator is anti-hermitian and projects the eigenvalues of -Dqv from the circle 



-16- 



onto the imaginary axis. At non-zero chemical potential Dp loses its anti-hermiticity and 
the eigenvalues come in three types: pairs of opposite real eigenvalues, pairs of complex 
conjugate imaginary eigenvalues and quartets of complex conjugate/opposite eigenvalues 
(A, —A, A*, —A*). To form the complex quartets at finite /i, typically two pairs of purely 
imaginary eigenvalues need to come together to allow for its creation; in our RMT discussion 
the former were called on-axis and the latter ojf-axis. 

The anti-unitarity symmetry of Dp follows from that of 75-Dw) Eq. (6.17), and the sign 
function thereof (and the facts that U commutes with 75 and that sign and 75 are their 
own inverses): 

75 sign*(75L>w) = sign(Z)w75)f^ = U'' sign{j5Dyy)'j5U , (6.20) 

D; = 2U^^ + ^ign(75Z^w)75 ^ ^ ^ [/t ^^M^M + 1 ^ 
^ 1 - sign(75L»w)75 75 sign(75Dw) - 1 

= -U^DpU. (6.21) 

Hence iDp shares the anti- unitary symmetry from the continuum, compare (6.21) to (6.7), 
and therefore GOE spacings are expected at vanishing Again a chemical potential 
destroys the hermiticity, but keeps the anti-unitary symmetry, exactly as for the continuum 
Dirac operator. Note that fj, enters the hermiticity breaking part of the operator Dp non- 
linear ly - unlike A and A in the random matrices. 

6.2 Numerical results 

The quenched lattice simulations were performed on an 8^ lattice with the Wilson gauge ac- 
tion, using the USQCD lattice QCD software Chroma and QDP-I— |- [47]. The software was 
adapted with an implementation of the overlap operator at non-zero chemical potential, 
based on the nested Krylov-Ritz approximation [48]. We generated 60000 quenched con- 
figurations at /3 = 2.2 using a heatbath/overrelaxation algorithm. On subsets of them we 
computed 0(30) of the lowest lying eigenvalues of the overlap operator (6.9) with ARPACK 
[49] for various values of the chemical potential and Wilson mass m„ = —1.4. The number 
of configurations used for the different chemical potentials are summarized in Table 1. An 
animated plot of the evolution of the eigenvalues of Dp with increasing chemical potential 
for a typical configuration can be found as ancillary material to this arXiv submission [50] . 

We measured the spacing distributions of the projected overlap operator Dp, defined in 
Eq. (6.19), for various values of the chemical potential As before, we distinguish between 
off-axis spacings and on-axis spacings with and without an interspaced complex pair. We 
only considered spacings between eigenvalues with an imaginary part in the spectral window 
(0.5,0.6). This ensures that the eigenvalues stay in the bulk of the spectrum (which of 





0.05 


0.10 


0.20 


0.30 


# configs 


60000 


30000 


20000 


20000 



Table 1. Values of the chemical potential and corresponding number of quenched configurations 
used to determine the spectral properties of the overlap operator. 
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0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 



s s s 




0.5 1 1.5 2 0.5 1 1.5 2 0.5 1 1.5 2 0.5 1 1.5 2 



s s s s 




0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 



s s s s 

Figure 5. Spacing distributions of the overlap operator with various real values of the chemical 
potential /x. Top: on-axis spacings of type 1, surmise given by GOE, Eq. (4.4). Middle: on- 
axis spacings of type 2, surmise given by Eq. (4.15). Bottom: off-axis spacings, surmise given by 
Eq. (4.24); the dashed curve is the GOE spacing for comparison. 

course is limited on a finite lattice) and have a roughly constant eigenvalue density which 
is necessary to apply the surmises, as argued in Sec. 5. As for the large random matrices, 
unfolding was only done for the on-axis spacings. 

The results are shown in Fig. 5. The numerical data are well matched by the surmises 
derived from small matrices, like for the large random matrices considered in Sec. 5. Again, 
the coupling parameter A is obtained by a fit for the on-axis spacings of type 2 and from 
the frequency p of imaginary spacings for the off- axis spacings. For the off- axis spacings at 
large fj, we observe that the discrepancy between the data and the surmise has the same 
tendency as for the large random matrices, cf. Fig. 4, lower right panel. This indicates that 
the RMT results for large matrices are able to describe the QCD results reasonably well, 
even at larger coupling. 

6.3 Imaginary chemical potential 

The introduction of an imaginary chemical potential does not change the anti-hermiticity 
of the continuum Dirac operator D nor of the lattice operator Dp. However, the operators 
iD and iDp no longer obey an anti-unitary symmetry. The random matrix ensemble ap- 
propriate for operators without anti-unitary symmetry is the GUE. Hence, with increasing 
imaginary chemical potential one expects a transition from GOE to GUE. Surmises for 
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Figure 6. Spacing distributions of the overlap operator with various imaginary values of the 
chemical potential ^ fitted by the surmise (6.22). For each ^, 5000 configurations have been analyzed 
in the spectral window with imaginary parts of the eigenvalues between 0.65 and 0.8. 



such mixed random matrices, again additive, Hqqe + A/Zgue, have been worked out in 
Refs. [12, 13], with the result 

Pgoe^gue(s) = Cse-^'^' erf ( ^) , (6.22) 



A 



where 



C = 2x/l + A2 D\ D = (r^-^ + arccotx) 



(6.23) 



In Fig. 6, these distributions are compared to those of the QCD Dirac operator Dp with 
imaginary chemical potential, were unfolding was again done by measuring the spacings in 
units of the local mean spacing (obtained by an ensemble average) . The coupling parameter 
was obtained by a fit with least square minimization. The surmise agrees very well with 
the numerical data. 



7 Summary 

We have shown that eigenvalues of real symmetric matrices are attracted under real anti- 
symmetric perturbations and are able to merge and move into the complex plain as complex 
conjugate pairs. For the emergent kinds of level spacings ~ on-axis with and without inter- 
spaced eigenvalue pair and off-axis - we have derived surmises from small random matrices. 
For the on-axis spacings without interspaced pair, the surmise is simply the GOE surmise, 
as in the unperturbed case, while for the other two we have obtained closed formulae, which 
depend on a coupling parameter A. 

These mixed- hermiticity surmises provide good approximations to the spacings of large 
random matrices in the regime of weak non-hermiticity and for the bulk of two-color QCD 
with small chemical potential, as expected from the anti-unitary symmetry and hermiticity 
properties of the Dirac operator. For on-axis spacings with an interspaced pair the surmise 
parameter A was obtained by a fit to the level spacings, while for off-axis spacings it was 
predicted by matching the frequency of those spacings between surmise and data. 

We have also measured the level spacings for two-color QCD with imaginary chemical 
potential and verified that they follow a mixed- symmetry surmise. 
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In contrast to the eigenvalue densities, the level spacings of asymmetric real random 
matrices have not been worked out analytically. However, for weak antisymmetric per- 
turbations the analytic surmises, which approximate the level spacings of large random 
matrices, can be used to describe physical systems obeying the relevant anti-unitary sym- 
metry properties. 
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